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I. INTRODUCTION 

As the lightest hadron, the Goldstone boson of spon- 
taneous chiral symmetry breaking and dominant parti- 
cle governing the nucleon-nucleon interaction, the pion 
occupies a special position in hadronic physics. Accord- 
ingly, tt-tt scattering has attracted considerable interest 
even though the cross section is not directly measurable. 
In particular, a variety of theories have been utilized to 
make scattering predictions: at low energies, current al- 
gebra/PCAC fjj, chiral perturbation theory §, |, |, § 
and bosonization models with pion fields ||; Quantum 
Hadrodynamic [QHD] meson exchange models and uni- 
tarized relativistic coupled channels Q| for intermediate 
energies; Regge theory at moderately higher energies; 
perturbative QCD at very high energies. These in turn 
have been confronted by experimental phase shift analy- 
ses H ||, [To| |ll| yielding reasonable agreement. However, 
there have been very few non-perturbative, relativistic 
covariant quark predictions which is the thrust of this 
paper. 

This work applies the quark Dyson-Schwinger [DSE] 
and meson Bethe-Salpeter [BSE] equations [[l2[ [b|, El) , 
in the rainbow-ladder approximation, to tt-tt scattering 
and extends our previous threshold analysis Jig] which 
reproduced Weinberg's low energy theorem |1| and the 
Adler zero. The approach is QCD based, renormaliz- 
able, relativistic, rigorously covariant, embodies cross- 
ing and chiral symmetry, and contains only two prede- 
termined parameters that provide a realistic, compre- 
hensive description of the light meson spectra and de- 
cays Jl|, H II . Our key finding is that the model 
correctly predicts all three isospin amplitudes measured 
at low energies and also reproduces the observed a and 
p resonances at higher energies. Our framework can thus 
assess phenomenological meson exchange models which 
require resonances as input and we find such models arc 
quite reasonable. 

This paper is organized into four sections. Section |0] 



formulates the scattering problem and addresses the main 
approximation: the rainbow truncation of the DSE for 
the quark propagator and the ladder truncation of the 
BSE for the pion as a quark- antiquark bound state. This 
section also details the quark annihilation and exchange 
diagrams and relates the amplitudes for these processes 
to the various isospin scattering amplitudes for total tt- 
tt isospin zero, one and two. Section III presents large 
scale, supercomputer calculations for forward and back- 
ward scattering and compares results to an effective me- 
son exchange model where the coupling constants and 
pole masses are a prerequisite. Finally, conclusions are 
summarized in section IV and supporting mathematical 
details are provided in an appendix. 



II. MESON SCATTERING IN THE 
DYSON-SCHWINGER APPROACH 

A. Dyson— Schwinger Equations 

The DSE for the renormalized quark propagator hav- 
ing four momenta p in Euclidean space is 



S(p) 1 = iZ 2 p' + Z 4 m,q(p) + 



Zi I g 2 D^{k)>^ ltl S{q)T a v {q,p), (1) 



where D^ u (k) is the dressed-gluon propagator, T"(q,p) 
the dressed-quark-gluon vertex with color index a = 
1...8, and k = p — q. The notation /^represents w^ki ■ 
The most general propagator solution of Eq. (Q) has the 
form S(p)^ 1 = i]/>A(p 2 ) + B(p 2 ) and is renormalized at 



spacelike p? according to A{p 2 ) — 1 and B(p 2 ) 



with m q (p) being the current quark mass. We use the 
Euclidean metric where {7 /t , 7„} = 2(5^, 7^=7^ and 

a ■ b = aibi = Ya=i 

Mesons are described by the Bethe-Salpeter amplitude 
[BSA], Th, which is a solution of the homogeneous BSE 
for q a q b bound states given by 
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K(p,q;P) 
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xs a ( q+ )r<i( q+ , q _)s b ( q -), (2) 

where a and b are flavor indices, p + = p + r/P and 
p_ = p — (1 — r?)P are the outgoing and incoming quark 
momenta respectively, and q ± is defined similarly. The 
kernel K is the renormalized, amputated qq scattering 
kernel that is irreducible with respect to a pair of qq 
lines. This equation has solutions at discrete values of 
P 2 = —m 2 H i where win is the meson mass. By then 
imposing the canonical normalization condition for qq 
bound states, Th is uniquely determined. Mesons with 
different spin, parity and C parity, such as pseudo-scalar, 
vector, etc., are characterized by different Dirac struc- 
tures. The most general decomposition for pseudoscalar 
bound states is p6[ 

r PS ( q+ , q -) = l5 [iE( q 2 , q -P ] r 1 ) + P F(q 2 ,q-P;r 1 ) 
+ ft G( q 2 , q ■P;rj)+ a MV P^ H( q 2 , q ■ P; r,)] , (3) 

where the invariant amplitudes E, F, G and H are 
Lorentz scalar functions of q 2 and q-P. For eigenstates of 
C parity, these amplitudes are appropriately odd or even 
in the C odd variable q-P. In the case of the h pion, for 
example, the amplitude G is odd in q ■ P, the others are 
even. Note also that these amplitudes explicitly depend 
on the momentum partitioning parameter r\. However, 
provided Poincare invariance is preserved, the resulting 
physical observables are -q independent |TJ|, [ij], [l?J] . 

B. Rainbow- ladder truncation 

We use the rainbow-ladder truncation for the system of 
Dyson-Schwinger and Bethe-Salpeter equations [DSBS]. 
In particular, the rainbow truncation of the quark DSE, 
Eq. (§, is 

z ig 2 D^(k)r:( q , P ) -> 5(fc 2 )D»7^, (4) 

where D^ c (k = p — q) is the free gluon propagator in 
Landau gauge, and Q(k 2 ) is an effective qq interaction 
that reduces to the perturbative QCD running coupling 
in the ultraviolet region. The corresponding ladder trun- 
cation of the BSE, Eq. (||), is 

K(p,q;P) - -g(k 2 )D f ™(k)^^- 7 „ , (5) 

where k = p — q . The two truncations combine to consis- 
tently produce vector and axial-vector vertices satisfying 
the respective Ward-Takahashi identities. In the axial 
case, this ensures that in the chiral limit the ground state 
pseudoscalar mesons are the massless Goldstone bosons 
associated with chiral symmetry breaking (l(| H- For 
vector mesons it yields a conserved electromagnetic cur- 
rent if the impulse approximation is used to calculate 
the electromagnetic form factor p9| . Furthermore, this 
truncation was found to be particularly suitable for the 
flavor octet pseudoscalar and vector mesons since the 
next-order corrections in a quark-gluon skeleton graph 
expansion significantly cancel p2] |. 



C. n-n scattering 

There arc different types of diagrams that contribute 
to 7T-7T scattering, or more general, meson-meson scat- 
tering, such as quark annihilation or quark exchange di- 
agrams. It is convenient to first discuss these different 
types of contributions in the impulse approximation, even 
though the impulse approximation is known to be insuf- 
ficient ||l5| . The more correct treatment is given in the 
next subsection. 

Figure [j] depicts the entire class of diagrams that con- 
stitute the impulse approximation for the scattering am- 
plitude. Each diagram can be associated with a specific 




FIG. 1: The six diagrams contributing to the impulse ampli- 
tude for n-n scattering; the incoming pions are labeled 1 and 
2, and the outgoing pions are 3 and 4. 

amplitude that is a function of the Mandelstam variables 
s = (Pi + P2) 2 , t = (Pi - P 3 ) 2 , and u= (Pi- F4) 2 , 
where Pi, P2 (P3, Pa) are the incoming (outgoing) pion 
momenta. The external pions are on-shell, yielding P 2 — 
—rri^ in the Euclidean metric. At threshold s = —4m 2 , 
t = u = 0, and the physical region has s < —Am 2 . Mo- 
mentum conservation requires s + 1 + u = —4m 2 7 so that 
the amplitudes, which we chose to express as functions 
of s, t, and u, only depend on two independent vari- 
ables. Alternatively the amplitudes could be represented 
as function of the pion center-of-momentum [cm.] en- 
ergy us = \J — s/2 and the cm. scattering angle 9 with 
cos((9) = (u — t)/(u + t). Then forward (backward) scat- 
tering, 9 = (9 = 180), corresponds to t = (u = 0). 

Because there are only two distinct microscopic scat- 
tering mechanisms, quark annihilation and quark ex- 
change, there are actually only two distinct amplitudes, 



D and E, respectively. Thus, all six diaj 


;rams in Fig. [l 


can be represented by these two 




diagram I 


D(s, t, u) 


(6) 


diagram II 


P(s, t, u) 


(7) 


diagram III 


D(s, u, t) 


(8) 


diagram IV 


D(s, t, u) 


(9) 


diagram V 


E(s,u,t) 


(10) 


diagram VI 


D(s, u, t) . 


(11) 
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These amplitudes can be calculated from integrals involv- 
ing the quark propagators and pion BSAs, e.g. 



D(s,t,u) = 2 N c / Tr S(k + Pi)T 7T (k + Pi, k)S(k) 

Jq 

x T„(k, k - P 2 )S(k - P 2 )f 1T (k -P 2 ,k + Pi- P 3 ) 
x S{k + Pi- P 3 )T„(k + P 1 -P 3 ,k + Pi)] . (12) 

Note that for tt-tt scattering, E is symmetric under inter- 
change of u and t (indicated hereafter by a semicolon), 
E(s;u,t) = E(s;t,u), whereas D is symmetric in s and 
t, D(s,t;u) — D(t,s;u). Furthermore, one can relate E 
to D: E(s;t,u) = D(u,t;s). However, we find it more 
convenient to retain and calculate both D and E inde- 
pendently since this affords a sensitive numerical check 
on our computer codes. 

The three isospin amplitudes, Tj, for tt-tt scattering 
are specific combinations of these amplitudes (see the 
appendix for further details). 

T (s,t,u) = 3(D(s,t;u) + D(s,u;t)) 

-±(E(s;t,u)+E(s;u,t)), (13) 

Ti(s,t,u) = 2(D(s,t;u)-D(s,u;t)) , (14) 
T 2 (s,t,u) = E(s;u,t) + E(s;t,u) . (15) 

Crossing symmetry and Bose statistics requires the even 
isospin amplitudes to be symmetric in exchange of u and t 
and only contain even partial waves. Similarly the isospin 
one amplitude must be antisymmetric in interchange of u 
and t and have only odd partial waves. Hence at thresh- 
old 4m^., 0, 0) = but the isospin zero and two am- 
plitudes Tq and T 2 are nonzero. 



D. Beyond the impulse approximation 

As documented in Ref . jjlj , the impulse approximation 
combined with the rainbow-ladder truncation is insuffi- 
cient to describe tt-tt scattering and it is necessary to 
include the ladder kernel as indicated in Fig. |[ Thus, 
for amplitude D we need to calculate and add an infinite 
series of s-channel ladder diagrams as well as an infinite 
series of t-channel ladder diagrams to the impulse term. 
Similarly for amplitude E, we must add and compute an 
infinite series of i-channel and u-channcl ladder diagrams. 

These infinite set of ladder diagrams can be calculated 



A 



FIG. 2: Diagrams needed to correctly describe tt-tt scattering 
in the DSBS approach in rainbow-ladder truncation. 



by solving an inhomogeneous BSE in ladder truncation 

,A 

F(p, P U P 3 ) = F (p, P^Pj) + / G(k 2 ) D*™(k) 

Jq 

where now q± — q ± (Pi + Pj) for i, j — 1,2 (s-channel 
ladders) or q± = q ± (Pi — Pj) for i — 1,2 and j = 
3,4 (t- and u-channel ladders). Here, F (p, Pi, Pj) is an 
inhomogeneous term of the type ST^; e.g. for the s- 
channel ladder contribution to configuration D we have 
i = 1, j = 2 and 

F (p,Pi,P 2 ) - T„(p + Q + ,p-Q_)S( P -Q-) 

xT 7r (p-Q_,p-Q + ) (17) 

where Q± = (Pi±P 2 )/2. This is similar to solving the in- 
homogeneous BSE for the quark-photon vertex, as done 
in Refs. |l8|, [l^]. The Dirac structure is indeed of equal 
complexity, in both cases there are eight independent 
terms. Thus, F can be written as 

F(p, p, Pj) = O n f n (p 2 ,P- PuP- Pj) (18) 



with 



1 (p,P i ,P j ) 
2 (p,P t ,Pj) 
3 (p,P t ,Pj) 

5 (p,P,Pj) 
6 (p,P,Pj) 
7 (p,P,Pj) 
Os{p,P,Pj) 



P 

Pi 

Pj 

PJ- Pi 
Pl-Pj 

Pi 1 ■ Pj 
PJ-PiJ-Pj 



(19) 
(20) 
(21) 
(22) 
(23) 
(24) 
(25) 
(26) 



Note that the /„ are functions of three independent vari- 
ables: p 2 , p-Pi andp-Pj, in contrast to the quark-photon 
vertex, which has only two independent variables. 



III. 



NUMERICAL RESULTS 



We consider only forward and backward scattering 
which simplifies the numerical analysis compared to scat- 
tering for arbitrary angles. However, since F is a func- 
tion of three independent variables for any given P\ and 
P 2 , there is a significant increase in the numerical ef- 
fort to solve Eq. ( |l6| ) when compared to the inhomoge- 
neous BSE for the quark-photon vertex. We discretize 
the three independent variables on a three-dimensional 
grid and solve by iteration, starting with Fq. If there 
are no singularities, the number of iterations for con- 
vergence is about 20 to 30. The actual calculations are 
done on a parallel processor supercomputer (IBM SP) 
and scale reasonably well, typically 2 to 3 CPU minutes 
on 32 processors per external Pi, Pj. Since there are two 



4 



sets of ladder diagrams for both configurations D and E, 
about 5 CPU hours per external momentum variable is 
required. The iterative procedure is significantly less effi- 
cient near a singularity, yet the entire parallel calculation 
still scales reasonable well. 



A. Model truncation 

We adopt the model framework and parameters that 
have been recently developed (TtJ which provide a good 
description of the masses and decay constants of the light 
pseudoscalar and vector mesons. The Ansatz for the ef- 
fective qq interaction is 

4tt 2 D k 2 



G(k 2 
k 2 



+ 



^ 2 lm T(k 2 ) 



I In 

- 2N f 



1 



k 2 /A 2 



QCD) 



(27) 



where 7m = f 2/(33 - 2N f ) and T(s) = (1 - exp ^)/s. 

The ultraviolet behavior is chosen to be the QCD run- 
ning coupling a(k 2 ); the ladder-rainbow truncation then 
generates the correct perturbative QCD structure for the 
DSBS system of equations. The first term implements 
the infrared strength in the region < k 2 < 1 GeV 2 
where the chiral condensate is fit (2j|. We use 
m t = 0.5 GeV, r = e 2 - f , N f = 4, A QCD = 0.234 GeV, 
and a renormalization scale p = 1 9 GeV which is in the 
perturbative domain jl^, |l7| . The remaining parame- 
ters, u = 0.4 GeV, D = 0.93 GeV 2 and the u/d degener- 
ate quark mass, are determined by fitting the chiral con- 
densate, and /„•. The predicted p meson mass and 
electroweak decay constant are also in good agreement 
with observation [jl7J , as can be seen from Table 0. Fur- 
ther, without parameter readjustment, the model agrees 
remarkably well with the most recent Jlab data [^4| for 
the pion charge form factor F 7T (Q 2 ) and the strong decay 



TABLE I: Calculated and measured meson masses and decay 
constant, adapted from Ref. |JL7[. 





experiment |25|] 


calculated 




(estimates) 


C fitted) 


„_ u — d 
m /i = lGeV 


5 - fO MeV 


5.5 MeV 


- m% 


(0.236 GeV) 3 


(0.241+) 3 




0.1385 GeV 


0.138 1 " GeV 


u 


0.924 GeV 


0.925+ GeV 


m p 


0.770 GeV 


0.742 GeV 


fp 


0.153 GeV 


0.146 GeV 



B. Results for isospin amplitudes near threshold 

We first elucidate the importance of the ladder con- 
tributions, reaffirming that the impulse approximation is 
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FIG. 3: Numerical results for D(s, t = 0; u) (top) and for 
D(s,t;u = 0) (bottom), using impulse approximation only 
(dotted), ladder contributions (dashed and dot-dashed), and 
total (solid), which is the sum of the impulse contribution and 
the two sets of ladder contributions. 



insufficient to describe 7r-7r scattering. In Fig. [| the im- 
pulse approximation (dotted curve) is compared to the 
contributions from the entire s-channel and t-channel lad- 
der diagrams for amplitude D(s,t;u) at extreme forward 
and backward scattering angles. The total amplitude 
(solid curve) is noticeably different than the impulse am- 
plitude, clearly indicating the ladder exchange diagrams 
are crucial. Also note that D(s, t; 0) is indeed symmetric 
in s and t as anticipated; within our numerical accuracy 
we also found that E(s; 0, u) = E(s; u, 0) = D(s, 0; u). 

Not only is the deficiency of the impulse approxima- 
tion clear, it is also essential that the full ladder of gluon 
exchanges be included, especially for kinematics near ir- 
7r resonances. This is displayed in more detail in Fig. |] 
where the effects from one, two, three, nine, and "in- 
finitely many" s-channel gluon exchanges are compared. 
To reach effective convergence requires up to about 20 
gluon exchanges away from the resonance region. How- 
ever, near intermediate bound state poles (p and a) con- 
vergence at the 1% level requires several hundred gluon 
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FIG. 4: Numerical results for contributions to D(s, t = 0; u) 
coming from the s-channel one gluon exchange diagram (dot- 
ted), one and two gluon exchange diagrams (dashed), one, 
two, and three gluon exchange diagrams (dot-dashed), up to 
nine gluon exchange diagrams (dot-dot-dashed), and a com- 
plete set of ladder diagrams (solid) . 



exchanges. 

Combining the amplitudes D and E with Eqs. (|l3|)- 
( |l5| ) provides the different isospin scattering amplitudes 
which are displayed in Fig. |^. For comparison we also 
plot the amplitudes from leading-order chiral perturba- 
tion theory |5] which are given in partial wave form 



Ti = 32vr^(2 L + 1)T 7 L (s)P L (cos ( 



(28) 



For low energies the S and P waves dominate and are 
given by 



T °(s) 



-2s- 



327T/2 

s + 2m\ 

327T/2 

— s — Ami 

967T/2 



(29) 
(30) 
(31) 



Again recall that in the Euclidean metric s — —Ami = 
-0.076 GeV 2 at threshold. 

Note the excellent agreement between tree-level chiral 
perturbation theory and our calculations at low energies. 
The deviation at higher energies represents contributions 
from 7r-7r resonances (scalar and vector meson bound 
states). These important physical effects are automat- 
ically included in the present approach and are discussed 
further in the next section. 

Finally, we have calculated the dimensionless S and P 
wave scattering lengths a\ in all three isospin channels 



= T°( S = -Ami) 



Ami 



Ami 



T?(s 



-Ami 



(32) 
(33) 



FIG. 5: The lines correspond to leading-order chiral pertur- 
bation theory for isospin zero (dashed), one (dot-dash), and 
two (solid) amplitudes near threshold. The symbols represent 
the DSBS amplitudes. 



We find ag = 0.17, a\ = 0.036, and a 2 , = -0.045. This 
agrees well with Weinberg's theorem which is embodied 
in the tree-level chiral perturbation results 



7< 

327T/2 

ml 

~2m 2 
32^/ 2 



0.156 



0.030 



= -0.044 



(34) 
(35) 
(36) 



A recent analysis pq] of the experimental data |27fl uti- 
lizes both two-loop chiral perturbation theory and a 
phenomenological description involving the Roy equa- 
tions fH and obtains ag = 0.220 ± 0.005, a\ = 0.0379 ± 
0.0005, and a% = 0.0444 ± 0.0010. Chiral perturba- 
tion theory is able to provide more accurate scattering 
lengths because they include higher order contributions 
from pion loops which is especially important for the 
isospin zero channel: one-loop chiral perturbation the- 
ory Q gives a[] = 0.200, compared to a[j = 0.156 at lead- 
ing order. In our present quark-based calculation pion 
loops are not included. However, they are clearly neces- 
sary for an accurate description of the data, in particular 
in the isospin-zero channel. 



C. Comparison with meson-exchange models 

It is enlightening to compare the DSBS approach to the 
QHD formalism which has an established phenomenolog- 
ical legacy. For example, one could formulate a meson- 
exchange model with both scalar and vector mesons 
starting from the effective Lagrangian 



Cj fL-x T Cstj T l~> p "F £h 



(37) 
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FIG. 6: Leading meson-exchange diagrams for n-n scatter- 
ing: (a) contact term, (b) s-channel, (c) i-channel, and (d) 
u-channel a and p exchange. 



Ant = 54tt ($n ■ $T^j + 9aTnr4'<T$TT ' 
+ 9/mTT<Pp ■ X ^p4>7T ■ 



(38) 



It entails meson masses and coupling constants which 
require phenomenological determination; it also requires 
a certain amount of fine-tuning in order to satisfy chi- 
ral constraints. The leading Feynman diagrams for 7r-7r 
scattering in such an approach are displayed in Fig. ^. 

In a meson-exchange model using both scalar and vec- 
tor mesons, the amplitude D(s, t; u) corresponds to a con- 
tact term plus meson exchange contributions in the s and 
t channels; the amplitude E(s; t, u) corresponds to the 
same contact term plus meson exchange contributions in 
the t and u channels. In order to calculate these me- 
son exchange contributions, let us first specify the meson 
coupling constants and gp^, which in impulse ap- 
proximation are given by 



9a 



4 Nr. / Tr 
'i 



s(k + p x )r w (fc + p u k)s(k)r n (k, k - p 2 ) 



x S(k-P 2 )T a (k-P 2 ,k + P 1 ) 



(39) 



4 N r . / Tr 



s(k + i\)r w (fc + p u k)S(k)r w (k, k-p 2 



xS(fe-p 2 )r;(fc-p 2 ,fc + Pi)] (40) 

for on-shell scalar and vector mesons, (Pi +P 2 ) 2 = —rn 2 , 
—m p respectively. 

Within the present model, we have m 2 = 0.449 GeV 2 , 
g„^ = 2.06 GeV, m 2 = 0.549 GeV 2 , and 5p7r7r = 
5.14 p3]. The corresponding width for the decay p° — > 



P7T7T 



9%* (m 2 p 



48 7T 

104 MeV 



4m 2 ) 2 



(41) 
(42) 



which is in reasonable agreement with the experimen- 
tal p width 151 MeV; the difference could very well be 
explained by pion loops which are not included in the 
present approach. The decay width for the a — > 7T7t is 



3 S 2 ™ \fi 



Ami 



2 16 7r ml 
= 172 MeV 



(43) 
(44) 



which is also quite reasonable for a broad a resonance 
(the factor | comes from summing over the charged and 
neutral pions). However, it is known that the scalar 
BSE receives significant corrections beyond ladder trun- 
cation [p2| , which could change this calculated decay 
width. In addition, pion loops should be incorporated 
self-consistently in the BSE approach for a more realistic 
calculation of the a and p widths. 

Next, consider for example the contribution of s- 
channel ladder diagrams to -D(s, t; it): scalar mesons con- 
tribute 



<7a™/4 



whereas vector mesons contribute 



9%* (« ~ *)/4 



(45) 



(46) 



near the bound state poles, and similarly for the t- and 
it-channel ladder diagrams. The momentum-dependent 
factor in the numerator of the vector meson exchange 
contributions comes from (Pi — P 2 ) ■ (P3 — P4) = u — t. 

A detailed, consistent evaluation of these diagrams 
generates the effective meson-exchange scattering ampli- 
tudes at tree level 



T (s,t,tt) 



1 



2 9<T ™ V s + ml ' t + ml ' u + m 2 



Ti(s,t,u) 
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9 pit IT I 7 l 9 ' l 9 
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(47) 



2 9plT7T V s + ml ' t + ml u + ml 
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(49) 



7 



We can now evaluate the accuracy of the meson ex- 
change model by comparison with the full DSBS am- 
plitude. A consistent assessment requires that the me- 
son parameters (i.e. the coupling constants and meson 
masses) must be provided by the DSBS model. The 
only free parameter in the meson exchange model is d n , 
which can be fitted such that the meson exchange model 
agrees with the microscopic calculation away from the 
resonances. A value of ~ — 4.8 provides reasonable 
agreement for both the isospin zero and two results over a 
large kinematic region. To reproduce Weinberg's results 
for the scattering lengths a[J and requires a fine-tuning 
between the meson masses and coupling constants in 
the meson exchange model. The meson exchange model 
with calculated m p , m a , g p7T 7r, and sw w from the DSBS 
model cannot reproduce the Weinberg limit: ctg = 0.156 
requires — —5.2, whereas a$ = —0.044 requires 

C^tt = —4.0. In contrast, the DSBS approach does re- 
produce the Weinberg limits correctly, while at the same 
time properly generating non-analytic effects from inter- 
mediate qq bound states such as p and a mesons. 



IV. CONCLUSION 

Summarizing, we have performed a consistent relativis- 
tic quark formulation of tt-tt scattering using the Dyson- 
Schwinger, Bethe-Salpeter framework. The DSBS ap- 
proach provides contact with QCD and naturally includes 
the important features of gauge and Poincare invariancc 
as well as both crossing and chiral symmetry. We have 
obtained numerical results which are in good agreement 
with the observed isospin scattering lengths and repro- 
duce Weinberg's theorem at threshold. Perhaps even 
more significant is the emergence of the a and p res- 
onances at the proper energy which reaffirms the QCD 
structure elements in the model. The latter permits a rig- 
orous assessment of QHD meson exchange models which 
simply insert resonances phenomcnologically. As we doc- 
umented in detail, a simple QHD treatment appears rea- 
sonable, especially near resonances. Future work will 
confront data at higher energies and also address other 
flavored systems such as n-K scattering. 
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FIG. 7: Isospin zero, one, and two ir-n amplitudes. The sym- 
bols are the DSBS calculations. The lines are meson exchange 
model results, using calculated meson masses and couplings 
from the DSBS approach. 



Fig. |?] summarizes our comparative study for all three 
isospin amplitudes at forward scattering. The gray band 
indicates the sensitivity on for —4.0 < < —5.2. 
Notice that the a resonance naturally appears in the 
DSBS S wave, isospin zero amplitude and agrees with 
the meson exchange amplitude quite accurately near the 
a mass region. Similarly the DSBS P wave, isospin one 
amplitude is reasonably well reproduced by the meson 
exchange result, although we cannot numerically quite 
reach the p pole. In between threshold and resonant 
regions the two approaches have some quantitative dif- 
ferences, however we conclude that the tree-level meson 
exchange approach is reasonable for simple predictions. 
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APPENDIX: ISOSPIN DECOMPOSITION 

In terms of the u, d flavor components, the single pion 
states are 



|i,i) = -M 



(A.l) 



|vr°) = |1,0) = \(uu-dd)/V2) (A.2) 
IO = = \du) . (A.3) 

The two-pion states then have isospin decomposition 



|7r ± 7r ± ) = |2,±2) 



I^tt ) = ^=(|2,±1)±|2,±1)) 



1 



V2 



(A.4) 
(A.5) 



k^) = -y2,0)±-yi,0) + ^=|0,0) (A.6) 



V6 



V2' 



k 7r ) 



(A.7) 



Using the latter equations, the physical 7r-7r scattering 
amplitudes can be related to the three isospin amplitudes 
by 



(tt^ISItt^) = T 2 



(A. 



6 T 2 + -T 1 + -T (A.9) 



8 



1 1 



1„ 1. 



(tt\ \S\it\ ) = 2 -T 2 + \ T{) 



(A.10) 
(A.ll) 
(A.12) 



Not all six diagrams in Fig. [I] contribute to each it- 
it scattering amplitude. In terms of the physical pion 
scattering amplitudes, we have 

(7r ± 7r ± |S'|7r ± 7r ± ) = E(s,t,u) + E(s,u,t) (A.13) 
(7r ± 7r =F |5'|7r ± 7r =F ) = 2 D(s,t,u) (A.14) 
(7r ± 7r°|,S'|7r ± 7r ) = D(a,t,u) - D(a,u,t) 



+ -(E(s,t,u) + E(s,u,t)) 



(A.15) 



(ti^ttFISItt ?!- ) = -D(s,t,u)- D(s,u,t) 



+ -(E{s,t,u) + E(s,u,t)) 



(A.16) 



(7T u 7r u |5|7r u 7r u ) = D(s,t,u) + D(s,u,t) 



+ -(E(s,t,u) + E(s,u,t)) 



(A.17) 



Thus we find for the isospin amplitudes 



T (s,t,u) = 3(D(s,t,u)+D(s,u,t)) 

-^(E(s,t,u)+E(s,u,t)) (A.18) 

Ti(s,t,u) = 2 (D(s, t, u) - D(s, u, t)) (A.19) 
T 2 (s,t,u) = E(s,t,u)+E(s,u,t) . (A.20) 

The quark exchange diagram, E(s]t,u), is symmetric in 
the last two arguments, t and u, whereas the quark anni- 
hilation diagram, D(s,t;u), is symmetric in the first two 
arguments, s and t. Furthermore, the quark exchange 
and quark annihilation diagram are related to each other 
via 



E(s]t,u) = D(u,t:s) . 



(A.21) 



Thus, there is actually only one independent amplitude, 
say D(s,t;u), where the symbol ; denotes the amplitude 
is symmetric in the first two arguments. In terms of this 
single amplitude, we have 



-D(u,t;s) (A.22) 
Ti(s,t,u) = 2 D(s,t;u) - 2 D(s,u;t) (A.23) 
T 2 (s,t,u) = 2 D(u,t;s) . (A.24) 

Finally, we recall the more conventional Weinberg am- 
plitudes A, B, and C defined by the standard scattering 
amplitude relation 

( 7 r 7 7r ,5 |S'|7r Q 7r' 3 ) = A{s, t, u)5 a p8^ 

+ B(s,t,u)S ai Sf3s + C^s, £, it)<S Q ,5<S/3 7 . (A. 25) 
Here the indices a, (3, 7 and 8 refer to the cartesian isospin 
projections (x, y, z) for the pion, which are related to the 
physical charge states by 



IT > = 7T 



(A.26) 
(A.27) 



An elementary isospin calculation immediately yields 

To(s,i, u) — 3 A(s,t, u) + B(s,t,u) 

+ C(s,t,u) (A.28) 

Ti(s,t,u) = B(s,t,u)- C{s,t,u) (A.29) 

T 2 {s,t,u) = B(s,t,u) + C(s,t,u) . (A.30) 



Then comparing with Eqs. ( |A~2^ )-( |A~2l| ) we can relate 
the three Weinberg amplitudes to the DSBS amplitude 
D 

A(s,t,u) = D(s,t;u)+D(s,u;t)-D(t,u;s)(A.3l) 
B(s,t,u) = D(s,t;u)-D(s,u;t) + D(t,u;s)(A.32) 
C(s,t,u) = -D(s,t;u) + D(s.u;t) + D(t,u;s) . 

(A.33) 

Since the Weinberg amplitudes satisfy 

B(s,t,u) = A(t,s,u) (A.34) 
C(s,t,u) = A(u,t,s) (A.35) 



we can also use Eq. (A. 31) to again derive B and C 



B(s, t, u) = D(s, t; u) + D(t, u; s) - D(s, u; t) (A.36) 
C(s, t, u) = D(u, t; s) + D(s, u; t) - D(s, t; u) (A.37) 



To(s,i,w) = 3 D(s, t; u) + 3 D(s, it; t) 



as an independent check. 
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